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We discuss several possible experimental signatures of the Bose-Einstein condensation (BEC) 
transition for an ultracold Bose gas in an inhomogeneous optical lattice. Based on the commonly 
used time-of-flight imaging technique, we show that the momentum-space density profile in the first 
Brillouin zone, supplemented by the visibility of interference patterns, provides valuable information 
about the system. In particular, by crossing the BEC transition temperature, the appearance 
of a clear bimodal structure sets a qualitative and universal signature of this phase transition. 
Furthermore, the momentum distribution can also be applied to extract the condensate fraction, 
which may serve as a promising thermometer in such a system. 
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I. INTRODUCTION 

There has been significant interest in ultracold atomic 
gases in optical lattices, partly stimulated by the possi- 
bility of simulating strongly correlated many-body sys- 
tems [ij. With extraordinary controllability, ultracold 
atomic gases in lattices provide a promising experimen- 
tal platform to help tackling many important problems in 
multidisciplinary fields. Among these topics, the emer- 
gence of condensation and superfluid order in an opti- 
cal lattice, and how the superfluid order transforms into 
other ordered states, is a problem which catches great 
attention over the past decade. With current technol- 
ogy, condensation and superfluidity are obtained for both 
Bose and Fermi gases in optical lattices, so that various 
phase transitions can be investigated [1, [3] ■ 

To experimentally investigate these phase transitions, 
present techniques with ultracold atomic gases heavily 
rely on detection based on the time-of-flight imaging, 
where the interference pattern and its visibility are sug- 
gested to be signatures of Bose condensation within an 
optical lattice [2]. Recent studies show that even for a 
thermal lattice gas above the BEC transition tempera- 
ture Tc, interference peaks with observable visibility are 
still present [1, 0, Q . This thermal visibility can be large 
for an ideal Bose gas in a homogeneous lattice, which 
could make the condensation signal ambiguous Q ■ How- 
ever, for practical systems with atomic interaction and 
an inhomogeneous global trap^he thermal visibility be- 
comes significantly smaller [5|, |6| , and the appearance of 
sharp interference peaks is still associated with the BEC 
transition. It is also suggested that the bimodal struc- 
ture of the atomic momentum distribution in the first 
Brillouin zone, combined with the interference peaks, 
provides an additional unambiguous signal for the Bose 
condensation in an optical lattice Q. A very recent ex- 
periment has used the onset of the bi-modal distribution 
and the associated condensate fraction to identify the 
supcrfluid-to-Mott-insulator transition point 

In this manuscript, we provide a detailed study of a 
Bose gas in a three-dimensional (3D) inhomogeneous op- 



tical lattice, both below and above the BEC transition 
temperature. We discuss several properties including the 
visibility, the width of the interference peak, and the mo- 
mentum distribution of the resulting interference pattern. 
The main results are as follows. First, all the quantities 
mentioned above can characterize the BEC transition for 
the experimental systems with interacting atoms in an 
inhomogeneous optical lattice. The large thermal visi- 
bility applies only to some particular parameters which 
are not directly responsible for current experiments. In 
the case when the thermal visibility is large, a substan- 
tial variation of the peak width or the appearance of a 
bimodal structure for the atomic momentum distribution 
may work as a better signature for the condensation tran- 
sition. Second, below the BEC transition temperature, 
the visibility and the peak width become insensitive to 
the system temperature, hence can not be applied as a 
practical thermometer. To fulfill this gap, the bimodal 
structure of the atomic momentum distribution gives a 
way to extract the condensate fraction through the bi- 
modal fitting. The resulting condensate fraction provides 
a sensitive indicator of the system temperature, hence 
may serve as a potential thermometer for this important 
system. 

The calculation techniques in this paper are similar 
to what we present in Ref. [6]. The remainder of this 
manuscript is organized as follows. In Sec. [TTl we first 
consider the situation of free bosons in an optical lattice 
within a global harmonic trap, and investigate the gen- 
eral behavior of the interference visibility and the atomic 
momentum distribution. In the absence of interaction, 
the problem is significantly simplified such that exact so- 
lutions are available. These exact solutions, on the one 
hand, are valuable for qualitative understanding of the 
system and its properties, and on the other hand, could 
be directly compared with experiments when the Fesh- 
bach resonance technique is applied to turn off the atomic 
interaction. After studying the free Bose gas, we then ex- 
tend our discussion in Sec. IIIII to the case of interacting 
bosons, where effects of the global harmonic trap and the 
interaction have to be taken into account together. In or- 



2 



der to deal with the trap, we adapt the local density ap- 
proximation (LDA), which works well when the interac- 
tion energy scale is significantly larger than the trapping 
energy scale (this condition is typically valid for current 
experiments) . Restricting our discussion to weakly inter- 
acting bosons that are away from the Mott region, we ap- 
ply the Hartree-Fock-Bogoliubov-Popov (HFBP) approx- 
imation to deal with the atomic interaction [11,, ,12. , ,13j . 
The HFBP method can provide a reliable description ex- 
cept for a narrow region around the BEC transition tem- 
perature [ll|, US • Taking into account the effect of the 
global trap, this questionable region only corresponds to 
a thin shell in three dimensions and its influence to the 
global properties is small. Therefore, we expect that the 
HFBP method can give reliable results for the atomic 
momentum distribution and the condensate fraction. 



II. IDEAL BOSE GAS IN AN 
INHOMOGENEOUS OPTICAL LATTICE 

In this section, we discuss an ideal Bose gas in an in- 
homogeneous optical lattice with a global harmonic trap. 
For completeness, we briefly review the formalism |6j] be- 
fore presenting various calculation results. We consider 
the atoms in a cubic lattice with an additional spher- 
ically symmetric harmonic trap [l5| . The Hamiltonian 
takes the form 



2m 



+ Vop{r) + V(r) 



(1) 



where ^ represents the bosonic field operator, m is the 
atomic mass, Vop(r) = Vb sin^(7rri/d) is the op- 

tical lattice potential with lattice spacing d, and V{r) = 
raitp'r^ 12 is the global harmonic trapping potential. In 
practice, the global harmonic trap ^(r) typically varies 
much slower than the optical lattice potential T4p(r), so 
the Hamiltonian can be separated into two parts with 
fast and slow variations, respectively. The fast-varying 
part can be diagonalized by introducing the expansion of 
bosonic field operators 



R 



w{r ~ R)aR, 



(2) 



where wiy) is the Wannicr function associated with the 
lattice potential Va^(y:\ or, is the annihilation operator 
on site R, and the summation is over all lattice sites. 
After transforming to the momentum space, the Fourier 
components of ^'(r), ti;(r), and or satisfy the following 
relation 



^'(k) = w(k)ak. 



(3) 



Representing the fast and slow varying components of 
H in terms of ak and or, respectively, the original Hamil- 
tonian Eq. ID) can be written as 



kelBZ 



R 



y(R)ai^aR, 



(4) 



where the summation over quasi-momentum k is re- 
stricted to the first Brillouin zone (IBZ). Here, we 
assume that the lattice depth Vq is strong enough 
such that the band gap is large and atoms are con- 
fined to the lowest band with dispersion relation ek — 
—'^tJ2i=x y z cos{kid). The tunneling rate t can be well 

estimated by i « (3.5/v^)V^o'^^'' exp(— 2vT^), where the 
recoil energy = H'^tt'^ / (2md'^) is used as the energy 
unit 

In principle, the resulting Hamiltonian Eq. Q can 
be directly diagonalized for arbitrary F(R). However, 
the numerical calculation is usually very heavy in three 
dimensions due to the presence of a large number of lat- 
tice sites. In this case, the diagonalization process can 
be significantly simplified by noticing that the indices k 
and R in Eq. ([4]) are reminiscent of the coordinate and 
the momentum variables in quantum mechanics. This 
observation allows us to write down a first quantization 
Hamiltonian corresponding to Eq. Q in the momentum 
space, where R is replaced by the momentum gradient 
— iVk @, H, B[l3|- The resulting effective Hamiltonian 
thus takes the form 



1 



(5) 



which represents free bosons with effective mass m* = 
h?/{m(jj'^) in a periodic potential ek with period |G| = 
211 /d along all three principal directions. Furthermore, 
since and ek are separable, this Hamiltonian can be 
reduced to three one-dimensional problems which require 
much less effort to solve. Notice that the properties of 
this effective Hamiltonian depend only on the ratio of 
h'^uj'^ / {tEn), which suggests that the variation of cu and 
Vq can be scaled to each other by fixing the dimension- 
less parameter fi^a;^/(V^^^'^e~^^^), where is used as 
the energy unit. In this section, we will keep Vq fixed 
(Vb = IOSr) and look at changes of the system properties 
under variation of the global trapping frequency lu. With 
different lattice barriers Vq, one can directly read out the 
result by simply re-scaling uj to keep h^uj'^ /{Vq^^c''^^^) 
fixed (a smaller barrier thus corresponds to a larger ef- 
fective trapping frequency). 

The quasi-momentum distribution (aj^ak) is then given 
by the square of the eigenstate wave functions 0n (k) of 
Hcs, where the expectation value is obtained by averag- 
ing over all eigenstates n with a Bose distribution factor 
g{En) — l/exp[/3(£'n — /i) — 1]. Here, En is the corre- 
sponding eigenenergy, is the chemical potential, and 
/3 — \/{kBT) is the inverse temperature. Taking into ac- 
count the presence of the Wannier function in Eq. (3), 
the atomic real momentum distribution is 



n(k) = (*t(k)vl/(k)) = |z«(k)p(atak) 
= \w{ktY.g{EnMnM'- 



(6) 



For a free gas, this momentum distribution remains un- 
changed during expansion, so that the signal from the 
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FIG. 1: Calculated column-integrated momentum density 
profile near the first Brillouin zone measured through the 
time-of-flight imaging, taking for two different trapping fre- 
quencies at various temperatures. Parameters chosen in these 
plots are close to those for a typical experiment of *^Rb, 
with lattice depth Vo — IOBr and a total number of par- 
ticles A'^ = 10^ . The global harmonic trapping frequencies are 
w = 27r X 10 Hz for (a)-(c), and 2-k x 60 Hz for (d)-(f). 



and increases by orders of magnitude when T crosses Tc- 
Since the total atom number satisfies ^ 1, the transi- 
tion is actually very sharp, and Tc is well defined by the 
requirement above. 

^From Fig. [l] one can see that for the case of a very 
weak global trap (w = 27r x 10 Hz), the interference peaks 
are indeed clearly visible even above the transition tem- 
perature. However, the BEC transition is still evident 
from the time-of-flight images as the interference peaks 
become much sharper when T gets below Tc- For the case 
with a stronger global trap (w = 27r x 60 Hz, which is close 
to the value in experiments), the interference peaks be- 
come blurred when T > Tc (although one can still read 
some pattern). Again, across the transition temperature 
Tc, the time-of-flight image, in particular the sharpness 
of the central peak, undergoes a dramatic change. 

These figures show that qualitatively a Bose conden- 
sation transition should be visible with the time-of-flight 
images. To have a more quantitative description, how- 
ever, it is desirable to have some single-value indicators 
which change sharply across the BEC transition so that 
one can characterize this phase transition by measuring 
the indicators. As possible candidates, we next discuss 
in detail two quantities, including the visibility of the in- 
terference pattern and the peak width associated with 
the atomic momentum distribution in Sec. Ill A1 and ITTB| 
respectively. While both of the two quantities can sig- 
nify the BEC transition in a reasonably strong global 
trap, the peak width becomes more accurate when the 
global trap gets weaker. After the condensation transi- 
tion, both of these two indicators become very insensitive 
to the variation of the system temperature, so they do 
not provide a good thermometer. Instead, we suggest to 
measure the condensate fraction from the bimodal fitting 
to the central interference peak, as discussed in Sec. Ill Cl 
The measured condensate fraction change continuously 
with the temperature, thus gives a good indicator for es- 
timation of temperature in this important system. 



A. Visibility of the interference pattern 



time-of-flight image taking along a crystallographic axis, 
say z direction, is just the columnar density 



[kx T^y) — 



i{k)dkz. 



(7) 



Using the technique sketched above, we show in Fig. [T] 
the calculation results for the column-integrated momen- 
tum distribution as record by the time-of-flight images, 
both below and above the BEC transition temperature 
Tc- In these plots, we choose parameters close to those 
for a typical experiment of *^Rb atoms, where the lat- 
tice depth Vo = lOii^R, and the total number of particles 
N = 10^. For this finite system, the transition temper- 
ature Tc is determined by requiring that the number of 
atoms in the ground state is of the order of 1 when T > Tc 



The visibility of the interference pattern has been in- 
troduced in Ref. [l3|. It is defined as the intensity con- 
trast of two characteristic points on the interference pat- 
tern [l3| 



(8) 



where 



and are (column-integrated) atomic intensi- 
ties at sites A and B, respectively, as shown in Fig. [ija). 
The point A represents the position of the secondary peak 
while B is along the circle of the secondary peaks where 
the intensity takes its minimum. The visibility defined in 
this way is clearly independent of the Wannicr function 
[the pre- factor in Eq. (3)]. 

The temperature dependence of the visibility v, espe- 
cially by crossing the BEC transition, is shown in Fig. [2] 
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FIG. 2: (color online) The visibility as a function of temper- 
ature around the BEG transition temperature Tc for various 
values of trapping frequencies. Other parameters are the same 
as those used in Fig. [T] with Vb = 10-Er and A'' = 10^. Remind 
that the variation of Vb is equivalent to that of u), with a fixed 
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for various values of trapping frequencies. For a very 
weak trap {lu — 2Tr x 10 Hz), the visibility is pretty high 
{v > 0.8) even with the temperature considerably larger 
than Tc- So in the limit of a vanishing u, this is con- 
sistent with the results in Ref. [1] for free bosons in a 
homogeneous lattice (without the global trap). However, 
for large trapping frequencies, as pointed out in Ref. 
the visibility becomes significantly smaller when T > Tc, 
leading to a more substantial drop across the BEC transi- 
tion. With a trapping frequency around w = 27r x 60 Hz, 
the visibility jumps should be pretty evident to observe, 
as shown in Fig. [21 However, the transition is stretched 
over a wide range of temperatures (the visibility begins 
to drop starting from a temperature significantly below 
Tc, see Fig. [2]), which may make the determination of the 
transition point from the visibility less accurate. Notice 
that with the scaling relation, a smaller barrier Vq corre- 
sponds to a larger effective trapping frequency uj. So with 
a shallower lattice, the change in the visibility across the 
BEC transition gets larger for free bosons. 

For the system temperature above Tc, there is no long 
range coherence in the atomic cloud, so the finite vis- 
ibility of the interference pattern is induced by residue 
short-range thermal correlations. To understand the dif- 
ferent behavior of the visibility, we calculate the short- 
range thermal correlation function around the trap center 
with different global trapping potentials. The real space 
correlation function is defined as 



C(R) = 



(9) 



with R = indicating the trap center. In Fig. [3l we 
show the correlation functions both below and above Tr 



FIG. 3: (color online) Gorrelation functions around the trap 
center indicating real space coherence below (for a superfiuid, 
SF) and above (for a thermal gas, TG) the BEC transition 
temperature, with trapping frequencies (a) lj = 27r x 10 Hz, 
(b) 2n X 60 Hz, and (c) 2n x 120 Hz. Gurves in each panel 
are taken at, from top to bottom, 0.9Tc(solid), 0.95Tc (solid), 
l.Orc(dotted), l.OSTc (dashed), and l.lTc (dashed), respec- 
tively. Parameters used in these plots are Vo = 10-Er and 
= 10^. For reference, the characteristic length of the sin- 
gle particle g round s tate wavefunction adjusted by effective 
mass (L = sJhJmFili) are (a) L = 23.1, (b) L = 9.4, and (c) 
L — 6.7, respectively. 



for different trapping frequencies. Notice that for weak 
trapping potentials [see, e.g., Fig. [Slja)], the correlation 
function extends to several lattice sites at temperatures 
above Tc, indicating the presence of thermal short-range 
coherence. As the trapping frequency increases, the cor- 
relation length for a thermal gas decreases as presented 
in Fig. mb) and [3fc) , which is consistent with the dis- 
appearance of interference peaks as shown in Fig. EJf). 



Momentum-space density profile and the peak 
width 



Up to now, we discuss the visibility characterizing 
the contrast of the interference pattern. In this subsec- 
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tion, we introduce another single- value quantity, the peak 
width, which characterizes the sharpness of the interfer- 
ence peak. We notice that while the visibility does not 
undergoes a sudden change across the BEC transition 
when the global trap is weak, the width of the central 
interference peak, however, always shrinks sharply when 
the condensation takes place. Therefore, the peak width 
is always a good indicator of the BEC transition inde- 
pendent of the strength of the global trap 

To introduce the peak width, first we look at the 
atomic momentum distribution, which gives more de- 
tailed information about the system. From Eq. (3), the 
momentum distributions in other Brillouin zones are sim- 
ply copies of the distribution in the first Brillouin zone 
weighted by the given Wannicr function, so it suffices to 
study the atomic momentum profile in the first Brillouin 
zone. In Fig. |4l we plot the column-integrated momen- 
tum distribution along one crystallographic axis (e.g., the 
X-axis) passing through the center of the first Brillouin 
zone. It is clear that while the density profile is a ther- 
mal distribution when T > T^, a bimodal structure starts 
to appear when a non-zero condensate fraction emerges 
in the lattice, characterized by a sharp peak at the cen- 
ter of the momentum space surrounded by a flat thermal 
distribution. The signal of this structural change is sig- 
nificant as soon as the system crosses the BEC transition 
[see Fig.Hb) and (e)]. 

In order to characterize the sharpness of the atomic 
momentum distribution in the first Brillouin zone, we 
introduce the peak width as a single-value parameter. 
For this purpose, we first define the middle value of 
n±{kx, ky) within the first Brillouin zone 



©/(27i)=10Hz 



a3/(27r)=60Hz 



,x 10 



_ 1 
^mid — 2 



max + mm 

(fex,fcy)eiBZ (fej;,fca)eiBZ 



n^{kx,ky). (10) 



In the simplest term, the peak width w is measured as 
the radius in the momentum space where the (column- 
integrated) atomic density nj_{kx, ky) first falls to this 
middle value rimid. In Fig. O we show the peak width 
as a function of temperatures for various trapping fre- 
quencies. Notice that the central peak width decreases 
monotonically with temperature, and most importantly, 
undergoes a sharp and substantial change by crossing the 
transition temperature. This distinctive feature is uni- 
versal for all trapping frequencies, hence can provide a 
clear criterion for the phase transition. However, as will 
be discussed later, the sharpness of the change of central 
peak width around Tc is guaranteed only for ideal Bose 
gases. In the presence of atomic interaction, the varia- 
tion of central peak width may be more flat and extended 
over a wider range of temperatures. 



C. The condensate fraction as a measure of 
temperature 

^From the discussion above, we notice that after the 
condensation transition, both the visibility and the peak 
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FIG. 4: (color online) Momentum space columnar density 
along the a;-axis in the first Brillouin zone. The left and 
right panels correspond to the cases of trapping frequencies 
uj = 2n X lOHz and 2n x 60Hz, respectively. Notice that while 
a thermal distribution is present above Tc [(c) and (f)], a clear 
bimodal structure starts to appear for temperatures slightly 
below Tc with only precentral condensate fractions [no ~ 0.14 
in (b) and ~ 0.09 in (e)]. Here, the solid curves represent the 
total momentum density profile, and the dashed curves repre- 
sent the momentum density profile of the normal component. 
Parameters used in these plots are Vo — 10-Er and A*" = 10^. 



width become almost flat to variation of temperature, 
as one can see from Figs. [2] and \5\ This means that 
it is hard to get any information about the temperature 
of the system from the measured values of visibility and 
peak width. Since temperature is one of the most im- 
portant quantity for the thermodynamical property of 
the system, it is desirable to have some experimentally 
measurable indicator which gives a good estimate of the 
temperature. The momentum density profile in princi- 
ple gives a lot of information, but it is not a single-value 
quantity which makes it hard to support a direct com- 
parison at different temperatures or for different systems. 
To overcome this drawback, we note that below Tc, the 
atomic momentum distribution in the first Brillouin zone 
always shows a bimodal structure, which actually gives a 
universal signal for the condensation transition. Further- 
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FIG. 5: (color online) The peak width (defined in the text) 
within the first Brillouin zone, taking across the transition 
temperature for various values of trapping frequencies. Notice 
that the sharp and substantial change around Tc occurs for all 
cases hence serves as a distinctive signature of BEG transition. 
Parameters used in these plots are Vb = IO-Er and A'^ — 10^. 



FIG. 6: (color online) The condensate fraction as a function of 
temperature below the BEG transition. By extracting no from 
the momentum distribution as shown in Fig.|3Ia-b) and (d-e), 
temperatures of such cases can be determined correspondingly 
(triangles) . Parameters used in this figure are — IO-Er and 
iV = 10^ 



more, from the measured momentum density profile, one 
can al"ways do a bimodal fitting to figure out the atomic 
fractions in the condensate and in the thermal parts, re- 
spectively. The measured condensate fraction thus can 
serve as a good estimate for the system temperature. In 
Fig. [6l we sho-w the calculated condensate fraction as a 
function of temperature for this system "with t"wo differ- 
ent global trapping frequencies. Notice that the conden- 
sate fraction no changes steadily and monotonically as 
the temperature T varies. ^From the relation no(T), one 
can estimate the temperature T through the experimen- 
tally measurable uq. Therefore, the condensate fraction 
no gives a single- value quantity which can serve as a cri- 
terion for the condensation transition (with no > 0); 
well as an indicator of the system temperature. 

III. INTERACTING BOSE GAS IN AN 
INHOMOGENEOUS OPTICAL LATTICE 

In this section, we discuss the more general and prac- 
tical case where the atoms in an inhomogeneous optical 
lattice has collisional interactions with each other. Many 
qualitative features discussed in the last section, how- 
ever, remain valid in the interacting case. The atomic 
interaction indeed brings up several different properties 
in a quantitative level. In this section, the emphasis of 
our discussion is on these differences. 

By taking the interaction into account, the time-of- 
flight images are expected to be modified in two major as- 
pects. First, the repulsive interaction between atoms will 
tend to broaden their spatial distribution in a trap, and 
hence narrow the corresponding momentum distribution. 
Second, during the time-of-flight expansion, the remnant 



atomic interaction transforms the interaction energy into 
the kinetic energy, in particular in the early stage of ex- 
pansion. As a consequence, the momentum distribution 
tends to be wider in the final images. The images also get 
a bit blurred due to the scattering of different interfer- 
ence peaks in the momentum space during the expansion. 
In this section, we will discuss both of these two points. 

For an interacting Bose gas in an inhomogeneous op- 
tical lattice, exact solutions as in the non-interacting 
case no longer exist. Instead, we use the local den- 
sity approximation (LDA) to treat the inhomogeneity in- 
duced by the global trap. When the interaction energy 
is much larger than the trap energy, each local region 
of the global trap behaves like a homogeneous system 
with interacting atoms in a pure optical lattice. The 
interaction in this local homogeneous lattice is then ana- 
lyzed under the Hartree-Fock-Bogoliubov-Popov (HFBP) 
scheme [III ^M- The validity of this approach is sup- 
ported by the following considerations. First, the LDA 
works well for a large number of atoms in a weak har- 
monic trap, which is the case for the parameter ranges 
considered below. Second, the HFBP method should be 
able to provide a reliable description of weakly interact- 
ing Bose systems, except for the region close to a phase 
transition (associated with the BEC or the Mott transi- 
tion). For an atomic gas in a global harmonic trap away 
from the Mott transition, such questionable region corre- 
sponds only to a thin shell in space, and its contribution 
to the global properties is far less significant. 

Under the LDA, we consider an interacting Bose sys- 
tem with a spatially dependent Hamiltonian 

H = ^{ei^-^i)alai,+j ^ a^qaLkflk'+qfl-k', (H) 

k k,k',q 
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where ek is the dispersion relation defined above, 
/i = /z(r) is the local chemical potential, and U = 
Uhg J \w{r)\^d'^r is the on-site interaction rate. For a 
typical experiment of ^^Rb, Uhg is related to the s- 
wave background scattering length as = 5.45 nm by 
Uhg = 4:Trh?as/m, and U takes an approximate form of 
U ~ 3.05Va^{as/d), with an energy unit of the recoil 
energy With the standard HFBP approach, we 

separate the bosonic operators into two parts: 



(12) 



where ■00 = {o-o) = (lo) represents the condensate com- 
ponent, and 6k is the fluctuation around it. After per- 
forming the substitution for at and into the original 
Hamiltonian Eq. terms that are cubic and quartic 

in (5^ and (5k will be present. These terms are reduced 
to quadratic forms under the HFBP by employing the 
Wick's theorem. As a result, we obtain a quadratic ef- 
fective Hamiltonian H^ff 



eff 



( eo - M + I no 

+ ^ [ck - + 2U (ntot - no)] (5^(5k 



k 

Uno 



Jk^-k + 44'5k 



(13) 



component, it should be emphasized that no leads to a 
delta function at zero momentum, which is an artificial 
result of LDA. In order to overcome this artifact, one 
needs to consider explicitly the broadening of the con- 
densate momentum distribution due to the presence of 
the harmonic trap. This can be done by the following 
procedure. First, we get the condensate fraction distri- 
bution no(r) over the trap. The condensate wave func- 
tion then can be well approximated by ipoil) ~ •\/'^o(r) 
under the Thomas- Fermi approximation [18]. The con- 
densate component of the quasi-momcntum distribution 
is thus given by the Fourier transform of the wave func- 
tion ■!/'o(r)- Second, by adding this condensate contribu- 
tion to that from the normal part, and multiplied by the 
Wannier function square |?«(k)p as in Eq. ([6]), we get 
the resulting momentum-space density distribution. 

Next, as in the case of an ideal Bose gas, we discuss 
several characteristics on the momentum-space density 
distribution. In Sec Hira and lllTBl we discuss the influ- 
ence of interaction on the atomic momentum distribution 
inside the trap, especially the associated visibility and the 
peak width, respectively. Then, we calculate in Sec. lIIICl 
the condensate fraction as a quantitative measure of the 
system temperature. Lastly, in Sec. IIII D| we analyze 
the interaction effect during the time-of-flight expansion, 
and conclude that this effect does not change much the 
characteristics discussed above. 



Here, eo = — 6t is the energy at the band bottom with t is 
the tunnelling rate defined above, tiq = ■(/'q is the per site 
density of the condensate fraction, and ritot is the total 
number of particles per site. In order to derive the ex- 
pression above, the saddle point condition is employed to 
make the coefficients of terms linear in b vanish, leading 
to the saddle point equation 



^(r) = eg — C/no + 2[/nt 



(14) 



This equation must be solved self-consistently with 
the number constraint ntot = — where = 
— (1//3) lnTr(e'''^'=") is the thermodynamical potential. 
This constraint leads to the number equation 



ntot = no + 



^ 2 

k^^O 



Ck 



to + Uno ^, f PE]^ 
- coth 



(15) 

where i?k = a/ (ck — eo + UnoY — U^n^ is the quasipar- 
ticle dispersion relation. 

By fixing the number density per site at the trap cen- 
ter, we can solve Eqs. (ITil) and (fT5)) self-consistently 
to obtain the chemical potential at the trap center fio- 
This result, together with the LDA relation of /i(r) = 
A*o — ^(r), allows us to calculate /i(r), and hence the 
condensate fraction no(r) and quasi-momentum distri- 
bution of the non-condensate part nk^o(r) at arbitrary 
position in the trap. The overall non-condensate quasi- 
momentum distribution, thus can be obtained by inte- 
grating over the whole global trap. For the condensate 



A. Visibility of the interference pattern 

We show in Fig. [7] the visibility of the interference pat- 
tern as a function of temperature across the BEC tran- 
sition. Since the presence of interaction sets another en- 
ergy scale, the system is not determined by the single pa- 
rameter of h'^Lo'^ /{tE^), as in the ideal gas case. Instead, 
we consider various combinations of lattice barriers and 
trapping frequencies. 

^From Fig. [Tl^a), we find that independent of the 
strength of the global trap, the thermal visibility with 
T >Tc remains small (with v about or below 0.4). This 
result is significantly different from the case of an ideal 
Bose gas, where the thermal visibility can be close to the 
unity for a very weak global trap. By crossing the BEC 
transition, the visibility clearly increases so that for an 
interacting Bose g high visibility signifies that the 
atomic cloud is in the condensate region. However, for a 
strong global trap (with a; ~ 27r x 120 Hz for instance), 
the variation of the visibility is not sharp at the transition 
point, but continues into a pretty wide region below T^. 
Thus, in this case, it becomes less accurate to use the 
jump of the visibility to determine the BEC transition 
temperature. 

Another feature we can read from Fig. [T^a) is the con- 
vergent behavior of visibility for a thermal gas with the 
same lattice potential but various values of trapping fre- 
quency. This behavior can be understood from the over- 
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FIG. 7: (color online) The visibility as a function of tem- 
perature around the BEG transition temperature Tc for (a) 
Vo = WEb. and (b) a; = 27r x 60 Hz. The number density per 
site at the trap center is set as ntot ~ 1, and the total number 
of particles in the trap is A'^ 10'*. 



all quasi-momcntum distribution 

Tik = J d^rni^{r) ^ An J drni^{r)r'^ , (16) 

where the integration is over the whole trap. Under the 
LDA, the spatial dependence of number density is only 
through the chemical potential nk(r) = rl,k(/^(r)), hence 
the integration can be rewritten as 



= 47r 



3/2 



/ £%(e)d\/e. 

■J un 



(17) 



For a thermal gas with a certain number density and 
chemical potential at the trap center, the function nk(e), 
and hence the integration over e in the equation above is 
fixed. Thus, all the trap can do is to re-scale the quasi- 
momentum distribution by a factor of w"^. For a given 
optical lattice characterized by a Wannier function, the 
momentum-space density profile of a thermal gas for var- 



ious values ofuj takes the identical shape, hence all signa- 
tures we can read from it must remain the same. Notice 
that since the LDA approach is reliable for all parameter 
ranges discussed here (with a; < 27r x 120 Hz), we con- 
clude that this result is the effect of a strong interaction 
compared to the trapping potential, which can not be 
smoothly connected to the non-interacting results. 

In Fig. [Zl^b) , we show the visibility for different lattice 
depths with a fixed strength of the global trap. With 
lower lattice depths, the thermal visibility increases. For 
Vo = 4£'r, for instance, the visibility varies almost lin- 
early with temperature near the condensation transition 
(for T from 0.5Tc to 1.2Tc), and a high thermal visibility 
remains (with v ^ 0.7 — 0.8) even when T crosses Tc- 
This fiat behavior hence makes it difficult to use visibil- 
ity to signify the condensate region and to identify the 
transition point for interacting atoms in a shallow lattice. 



B. Momentum-space density profile and the peak 
width 



After discussing the visibility of interference peaks, we 
next focus on the momentum distribution in the first Bril- 
louin zone. In Fig. [51 we show the column- integrated mo- 
mentum density profile along one of the crystallographic 
axis (say, the a;-axis). Similar to the case of an ideal 
Bose g clear bimodal structure appears for tempera- 
tures below Tc even when the condensate fraction is still 
small [see Fig. |5l^b) with uq ~ 9%]. Therefore, the mo- 
mentum distribution and its bimodal structure sets an 
unambiguous criterion for the BEC transition, especially 
when supplemented with the interference pattern from 
the lattice structure. 

To characterize the sharpness of the central peak, we 
still use the peak width defined above as a single-value 
parameter. As shown in Fig. [51 the peak width signifi- 
cantly reduces across the BEC transition, so a sharp peak 
with a small width sets a clear indicator that the system 
is in the condensate region. This conclusion is qualita- 
tively consistent with the case of an ideal Bose gas as 
discussed in Sec. IIIBI For an ideal Bose gas, the peak 
width always has a sharp and large jump at the BEC 
transition point. For an interacting Bose gas, this jump 
becomes less sharp under certain circumstances. From 
Fig. [9l we notice that the jump of the peak width re- 
mains sharp across Tc for weak global traps or for lower 
lattice depths. However, in a stronger global trap with a 
higher lattice depth, the decrease of the peak width takes 
place over a wider range of temperatures, which makes 
it less accurate to locate the transition point using the 
peak width as an indicator. If one compares Figs. [3 and 
[51 it is interesting to note for that for a weak optical lat- 
tice, even when the visibility becomes too flat to show 
a phase transition, the variation of the peak width re- 
mains significant and sharp to serve as an indicator of 
the condensation transition. 
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FIG. 8: (color online) Momentum space columnar density 
along the x-axis in the first Brillouin zone, where clear bi- 
modal structures appear for temperatures below Tc. Here, 
the solid curves are the total momentum density profile, and 
the dashed curves are the momentum density profile of the 
normal component. The trapping frequency used in these 
plots is = 27r X 60 Hz, and the number density per site is 
unity at the trap center. 
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C. The condensate fraction 

As in the case of an ideal Bose gas, when one moves 
into the condensate region, the visibihty and the peak 
width become insensitive to the variation of tempera- 
ture. Instead, we need to use the condensate fraction 
as an indicator of the temperature in the condensate re- 
gion. The condensate fraction can be measured similarly 
through the bimodal fitting to the atomic momentum 
distribution. In Fig. [TOl we show the total condensate 
fraction as a function of the temperature under a couple 
of different lattice barriers. First, the condensate frac- 
tion change monotonically and sensitively with tempera- 
ture, so it provides a potentially good thermometer. Sec- 
ond, it is interesting to note that the condensate fraction 
does not approach the unity even when the temperature 
tends to zero. This result is significantly different from 
the case of an ideal Bose gas, where at zero tempera- 
ture it is always a pure condensate. This discrepancy 
is due to quantum depletion of the condensate at zero 
temperature, which is always present for an interacting 
gas. In the zero temperature limit, the HFBP approx- 
imation used here reduces to the Bogoliubov approach, 
so it naturally takes into account the contribution from 
quantum depletion. In Fig. Illl we show the in-trap zero 
temperature quantum depletion fraction for various lat- 
tice barrier depths. Notice that by applying a higher lat- 



FIG. 9: (color online) The peak width within the first Bril- 
louin zone, taking across the transition temperature for (a) 
Vo = IO-Br and (b) oj = 27r x 60 Hz. The number density per 
site is unity at the trap center. 



tice barrier Vb , the quantum depletion fraction gets more 
significant. This observation justifies the zero tempera- 
ture results of condensate fraction as shown in Fig. [TUl 
and is consistent with the trend one should expect. In 
fact, in the case with no optical lattice, the condensate 
fraction should approach unity at zero temperature for a 
weakly interacting dilute gas (with a small gas parame- 
ter). In the opposite limit, when the lattice barrier tends 
to the critical value for the Mott transition, the conden- 
sate fraction should deplete to zero. When the quantum 
depletion is dramatic, the HFBP is no longer a good ap- 
proximation. However, for parameter ranges discussed 
above, even with Vb — IZE'r, the condensate fraction 
(~ 70%) still dominates at zero temperature, assuring 
the validity of the HFBP approximation used through- 
out this section. 

Another interesting feature of the condensate fraction 
plot Fig. [in] is the crossing behavior of the two curves for 
different Vq at finite temperature. For a homogeneous gas 
with a certain number density, one would expect a more 
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FIG. 10: (color online) The condensate fraction as a function 
of temperature below the BEG transition, for an interacting 
Bose gas in an optical lattice with a; = 27r x 60 Hz and different 
lattice barriers Vo- The number density per site is unity at the 
trap center. By comparing with the momentum distribution 
as shown in Fig. [S] temperatures of such a system can be 
determined correspondingly (triangles). 




FIG. 11: The in-trap quantum depletion fraction as as a 
function of the lattice barrier Vo- The trapping frequency is 
= 27r X 60 Hz. The number density per site at the trap 
center is unity. 



severe depletion (both quantum and thermal) from the 
condensate for the higher Vq case, since both temperature 
and interaction would play relatively more significant role 
when the transfer integral t decreases. However, this sim- 
ple trend becomes more complicated in a trapped case. 
In fact, due to the presence of a global trap, thermal 
wings must emerge at the edge at any finite tempera- 
ture. As one calculate the in-trap condensate fraction, 
an integration over the whole trap should take the ther- 
mal wings into account, and the result will inevitably 
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FIG. 12: (color online) The number and condensate density 
distribution along the radial direction away from the trap 
center. Two different temperatures are considered as in (a) 
T = O.lTc, and (b) T = O.ST^. The solid curves correspond 
to the overall density profiles and the dashed curves to the 
condensate portions. The trapping frequency used here is 
liJ = 27r X 60 Hz. The total number of particles A^tot and 
the number of condensate particles A^o are respectively (a) 
Vo = 6SR(red): iVtot = 0.54 x 10*, iVo = 0.50 x 10*; 'Vo = 
10£;R(black): iVtot = 1.3 x 10*, A^o = 1.0 x 10*; (b) Vo = 
6£;R(red): iVtot = 2.5 x 10*, iVo = 2.8 x 10^; Vo = WEr 
(black): iVtot = 2.3 x 10*, No = 6.7 x 10^ 



depend on how much the thermal part contribute to the 
total number. As an example, we show the in-trap num- 
ber and condensate density profiles in Fig. [12] with fixed 
total number density at the trap center. At a lower tem- 
perature T = O.lTc, the thermal contribution is negligi- 
ble so it is apparent that a higher depletion is observed 
in higher barriers, as shown in Fig. I12r a). At a higher 
temperature T = O.STc, the thermal part becomes sig- 
nificant and one observes a more extended thermal tail 
in the lower barrier case with Vq = 6i?R, as shown in 
Fig- fTSl b) . This greater thermal contribution to the to- 
tal number makes the condensate part relatively smaller 
and leads to a smaller condensate fraction. Therefore, 
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the crossing behavior of no as shown in Fig. [10] is a con- 
sequence of the inhomogeneous distribution of different 
phases with various number density and chemical poten- 
tial, which are results from the inhomogeneity of the trap. 



D. Interaction effects during the time-of-fiight 
expansion 

Up to now, we have calculated the atomic momentum 
distribution inside the trap and have neglected modifica- 
tion of this distribution caused by the atomic interaction 
during the time-of-flight expansion. During the expan- 
sion of the atomic cloud, the initial interaction energy is 
transferred to the kinetic energy, which tends to broaden 
the momentum distribution [l9[. This influence is most 
evident for the condensate part, due to its high number 
density and narrow momentum distribution at the be- 
ginning. The influence on the momentum distribution 
of the thermal and the non-condensate part is negligi- 
ble, as confirmed by experiments [5], since this part has 
a broad initial momentum distribution already, and the 
weak atomic collisions are unlikely to cause any signifi- 
cant modification. In the following, we only discuss in- 
fluence of the atomic collisions on the condensate part of 
the momentum distribution. 

After turnoff of the optical lattice and the confining 
potential, the atomic cloud undergoes a free expansion 
in space, and the expansion of the condensate part can 
be well described by a time dependent Gross-Pitacvskii 
(TDGP) equation 

ihdt^iv, t)^(^^^ + U\^{r, t)\'^ ^P(r, t), (18) 

where the initial condition ^'(r, i = 0) {t is the expansion 
time) is given by the equilibrium condensate wave func- 
tion inside the optical lattice and global trap, which has 
been calculated with the method detailed in the above. 
Note that only for the condensate part we use the TDGP 
to evolve its momentum distribution. 

It is easier to understand the consequence of this evolu- 
tion by looking at the TDGP equation in the momentum 
space. The Fourier transform of Eq. gives 

Zm 

+ *-k+qW*-k'+q(0*k'(0- (19) 

k',q 

Clearly, without atomic collisions (the U term), the mo- 
mentum distribution |^'k(^)P remains unchanged. The 
atomic collisions transfer a pair of atoms from momenta 
(k', — k' -l-q) to (k, — k + q), which modulates the overall 
momentum distribution. 

To understand the consequence of the collision induced 
momentum transfer, we note that the initial condensate 
momentum distribution has a number of peaks, one in 
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FIG. 13: Illustration of the scattering between atoms from 
different peaks (the central and the secondary). This process 
leads to blurring of the interference peaks with a characteristic 
blurring pattern. 



each Brillouin zone. As the central peak in the first Bril- 
louin zone (with k close to zero) is the highest one, the 
scattering of pairs of atoms around the first peak satis- 
fying the momentum conservation has the largest contri- 
bution to the collision effect. Around the central peak, 
the wave function ^'k(i) has approximate spherical sym- 
metry, and the evolution of the momentum distribution 
around the central peak by the TDGP equation has been 
calculated and shown in Ref. The condensate peak 
gets somewhat lower and broader, however, its width is 
typically still significantly less than the width of the ther- 
mal cloud, and a bimodal structure of the momentum 
distribution remains clearly visible. Thus, the collision 
during the time-of-flight expansion has some quantitative 
influence on the peak width we calculate before, but the 
effect is not large and it should not change all the qualita- 
tive discussions in the last sections. In particular, as the 
bimodal structure remains clearly observable, we do not 
expect that the condensate fraction measured through 
the bimodal fitting to the momentum distribution has 
any significant change by the collision effect during the 
expansion. 

In the next order, the atom in the central peak can col- 
lide with another atom in the secondary peak, and scatter 
to some other directions in the momentum space. The 
momentum difference between the central peak and the 
secondary peak is given by G (G is related to the lattice 
constant d through |G| = 27r/(i), so the kinetic energy 
difference between them corresponds to a large energy 
scale fi.^|Gp/(2m), which is typically larger than the in- 
teraction strength during the expansion [the latter can 
be estimated by Un (i), where n (t) is the instantaneous 
atomic density at the collision]. Therefore, to be effec- 
tive, the collisions need to satisfy the momentum con- 
servation as well as an approximate energy conservation 
in the momentum space. As a consequence, for atoms 
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with incoming momenta around and G (corresponding 
to the central and the secondary peaks, respectively), the 
outgoing atoms are centered around a spherical surface 
in the momentum space as shown in Fig. [12] [with mo- 
menta (k',-k' + G), where k'^ + (-k' + G)^ « |G|2]. 
The sphere has origin at k = G/2 and a radius of |G|/2. 
This collision effect causes some blurring of the original 
peaks. Since the scattered atoms are dominantly around 
a sphere in the momentum space, this blurring scatter- 
ing gives some characteristic momentum distribution pat- 
tern. Experimentally, by looking at such a pattern, one 
may measure and constraint the magnitude of collision 
effects during the time-of-flight expansion. 

At high orders, there could be also scattering between 
different secondary peaks as well as scattering between 
the central and even higher order peaks. Although these 
scattering can change some of quantities we calculate be- 
fore, we expect all these modifications are pretty small. 

IV. CONCLUSION 

In summary, we have discussed in this manuscript ideal 
as well as interacting Bose gases in an inhoniogcneous op- 
tical lattice within a global harmonic trap. By explicitly 
calculating the momentum distribution, we have stud- 
ied several possible signatures of the BEC transition in 
an lattice based on the common detection technique of 
time-of-flight imaging. For parameters of relevance to the 



current experiments, a large visibility, a substantial de- 
cline of the peak width, and the appearance of a bimodal 
structure for the central peak, can all be used as signals of 
the condensation transition as one decreases the temper- 
ature. For some other parameters, the thermal visibility 
could be significant, and in such a case the latter two 
criteria will work better. In particular, the appearance 
of a bimodal structure for the momentum distribution 
is a robust signal associated with the condensation tran- 
sition in both free space and lattice (in the lattice case, 
the interference peaks give further information about the 
underlying lattice structure). 

After the condensation transition, both the visibility 
and the peak width become insensitive to the variation 
of temperature, so they can not serve as a practical ther- 
mometer. Instead, one may measure the condensate frac- 
tion by a bimodal fitting to the atomic momentum dis- 
tribution. The condensate fraction changes steadily with 
temperature, and may work as a good experimental indi- 
cator of the system temperature by comparing with the 
results from theoretical calculations. 
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